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Abstract 

We show some improved mapping properties of the Time Domain Electric Field 
Integral Equation and of its Galerkin semidiscretization in space. We relate the weak 
distributional framework with a stronger class of solutions using a group of strongly 
continuous operators. The stability and error estimates we derive are sharper than 
those in the literature. 
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1 Introduction 

The electric field integral equation (EFIE) has received much attention as a representation 
of electric fields that is conducive to efficient discretization schemes. The EFIE and its 
variants are a competitive alternative to the finite element discretization of the Maxwell 
equations due to two main advantages: (a) reducing problem dimension by one and, es¬ 
pecially, (b) handling unbounded domains in a natural way. However, concerns over its 
numerical stability hamper the method’s popularity in practical applications. Several pa¬ 
pers, for example [131 EZl UHl E], address this issue with various techniques, yet most of 
them focus on the frequency domain analysis. Important as it is, the frequency domain 
equation is ill-suited to deal with broad band waves, which can be adequately treated 
in the time domain. Rigorous numerical study of the Time Domain Electric Field Inte¬ 
gral Equation (TDEFIE) is, to our knowledge, scarce in the literature. Early attempts 
[2311I1IS] all develop the kind of frequency domain analysis debuting in EE] and then 
obtain the time domain estimates by inverse Laplace transforms or Plancherel identities. 
A lone contribution using time-domain techniques for the analysis —not dealing with 
discretization, and limited to smooth domains— is given in es- 

Pure time domain analysis has been shown to outperform the double-back-through- 
Laplace-domain approach in several situations. Time domain analysis in this context 
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originated in [23] as a tool to analyze long term stability of several boundary integral 
formulations in acoustics, and was developed in [12] to provide improved bounds for the 
retarded layer potentials and integral operators of transient acoustics. The same approach 
was further developed and refined for the direct integral equation formulation of transient 
scattering by a sound-soft obstacle [1], a boundary integral formulation for transmission 
problems in acoustics [20] , and indirect formulations for Dirichlet and Neumann problems 
for acoustics [21]. What distinguishes this paper from the previous time-domain analysis 
is the transformation of the abstract second order differential equation associated to a 
dissipative operator to a system of first order equations (in time as well as in space). 
This simplihes the analysis with respect to [nmiEQ] by avoiding the introduction of a 
cut-off boundary that was required to £t the problem in the correct functional framework. 
We note here that the techniques of [21] are not easy to adapt to the TDEFIE, due to 
non-minor complications in the associated Sobolev spaces. 

The results in this paper improve the Laplace domain estimates in two ways: (a) less 
regularity is needed of the input data, i.e., the mapping properties are improved; (b) the 
constants in our estimates are independent of time, leading to a reassuring conclusion that 
the solution will not blow up as long as the input data is compactly supported. We also 
carry out the analysis meticulously to reveal the constant’s dependence on the velocity of 
the wave, which might eventually shed some light on the low frequency breakdown of the 
EFIE. Remarkably, these results are valid irrelevant of the scatterer’s shape regularity, 
be it smooth or polyhedral, all thanks to the foundational work 0I7|. They naturally 
incorporate as a special case the mapping properties of the boundary integral operators 
and layer potentials. We envision our results to be a unifying and instrumental step in 
carrying out the analysis of the full discretization for different time semidiscretization 
strategies like finite differences, space-time Galerkin or convolution quadrature. It is 
also surprising (although this is not new [H [9]) to note that the TDEFIE is amenable 
to a general Galerkin discretization-in-space, and no discrete Hodge decompositions are 
needed, as opposed to the requirements of the frequency domain EFIE [T2] . 

The paper is structured as follows. In Section 2, we write the time-dependent Maxwell 
equations into the TDEFIE and state the main results, concerning the regularity of the 
solution as a function of the data, before and after Galerkin semidiscretization-in-space, 
error estimates for Galerkin semidiscretization, and mapping properties of the forward 
operator. In Section 3, we set up a mathematically rigorous background to understand 
the TDEFIE equation in the distributional sense. In Section 4, we apply Gq— semigroup 
theory techniques to prove the main results, which are then compared with the existing 
results in the literature in Section 5. 

Foreword on notation and background. Standard Sobolev space notations like 
L^{0),H’^{0) for boundary or domain O are assumed throughout the paper. Bold¬ 
face notation, such as L^((T) and H™'((T), is used for vector valued functions with each 
component in the corresponding scalar function space. Given an open set O, || • Ho de¬ 
notes the norm of both L‘^{0) and L^(G), while (•, ■)o is the associated inner product. We 
will write Y) to denote the space of bounded linear operators between two Banach 
spaces X and Y. For some very basic background on causal vector-valued distributions of 
a single variable and their Laplace transforms, we refer to m Ghapters 2 and 3], which 
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reduces the scope of a vaste theory that can be explored in [26] or [TT]. The space of of 
fc-times continuously differentiable functions from an interval J C M to a Banach space 
X is denoted 


2 Main results 


Let C be a bounded Lipschitz open set with boundary T, let f2_|_ := \ be its 

exterior, and let jz : T —)■ be the outward pointing unit normal vector held on T. Our 
goal is the study of the properties of a boundary integral formulation for the Maxwell 
equations in vacuum, written in terms of the electric held, in the exterior of a perfectly 
conducting scatterer occupying f2_. In PDE notation, the problem can be written as the 
search for E ; 0_|_ x [0, cxo) —)■ (the scattered electric held) such that 


'^E + VxVxE = 0 

in 12+ X [0, cxo). 

(2.1a) 

(E + E‘“) xiy = 0 

in T X [0, oo). 

(2.1b) 

o 

o 

in 12+, 

(2.1c) 

atE(-,0) = 0 

in 12+. 

(2.1d) 


Here is a given incident held and the curl operator is taken on the space variables. The 
scattered held is assumed to be radiating, i.e., E(-, t) has a bounded spatial support for all 
t. We will use a single layer representation of E in terms of an unknown J : T x R —)■ 
satisfying 

J(-,f) = 0 Vf < 0, J(-,f)x^' = 0 Vt, (2.2) 

that is, J is a causal tangential vector held on T. The single layer potential ansatz is given 
by the formula [IHl 2.2.3] 


E(x,t) = -c ^ J 


+cV 


47r|x — y| 
1 

47r|x — y| 


■atJ(y,t-c-i|x-y|)dr(y) 


(2.3) 


divrJ(y, T)dT dr(y), (x, i) e x [0, oo), 


where divp is the tangential divergence operator, and V is the gradient in the x variable. 
The tangential component of the Maxwell single layer potential on points x G T is 


(Vc*J)(x,t) : = 


-c X 


47r|x 




-ii 


y|)dr(y) X 




+c Vr 


47r|x — y| 


divrJ(y,r)dr 


j dr(y) j , 


where Vr is the tangential gradient. We emphasize that so far our presentation is merely 
formal and the convolution symbol in the dehnition of Vc * J is just notation. The single 
layer potential will be denoted E = «Sc * J and it is dehned on both sides of the boundary. 
The tangential component of the single layer potential is continuous accross the boundary, 
namely 

(E~ - E+) X V = 0, 
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where the ± superscripts denote limits from The time derivative of the density can 
be recovered from tangential values of the single layer potential with the formula 

dti = c (V X E'^ — V X E“) X u. 

More details on the meaning of tangential trace operators will be given below. The only 
reading of the incident held on the boundary of the conductor is given by its tangential 
component: 

/3(x, f) := — ^'(x) X E™‘^(x, t)|r x i>'(x). 

The Time Domain Electric Field Integral Equation (TDEFIE) determines J by imposing 
the boundary condition 

(Vc * J)(x,f) =/3(x,t), (x,t) e r X [0, oo). (2.4) 

The solvability analysis for this equation will be carried out at the same time as the 
analysis of a semidiscrete-in-space version of it. Galerkin semidiscretization starts with a 
hnite dimensional space Xh containing elements satisfying: 

G L2(r), X 1/ = 0, divr/x'^ G ^^(r). 

For instance, if T is a polyhedron, the Raviart-Thomas or Rao-Wilton-Glisson |2T] ele¬ 
ments dehned on a triangulation of T can be used as a discrete space. The semidiscrete 
Galerkin equations look for : T x M —)■ satisfying 

j'*(-,f) = 0 f<0, t>0, (2.5) 


and 


(Vc * J'')(x,f) ■ /x^(x)dr(x) 


/3(x, f) ■ /x^(x, t)dr(x) G X/i, Vf. (2.6) 


Jr Jr 

The approximated electric held is the result of inputting the density in the single layer 
potential expression E^ = «Sc * J^. 


Weak tangential traces. In order to state the main result of this paper we need some 
Sobolev space notation. For an introduction to Sobolev spaces related to the Maxwell 
equation we refer to [H]. Trace theorems, and a full characterization of the range of 
the trace operators, are studied in [5], |6] for polyhedra, and in [7] for general Lipschitz 
domains. We briehy recall some dehnitions and results. In the space 

H(curl,D+) := {n G L2(D+) : V x n G L2(D+)}, 

endowed with its natural norm 

ll^ll™rl,0+ •= ll^llo+ + 11^ ^ ^llo+) 

we can dehne a tangential trace operator that extends = u x ix. This operator is 
bounded and surjective from H(curl, 12+) to H~^/^(divr, T), which, roughly speaking, is 
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the space of H tangential vector fields whose surface divergence is in H 

The space H~^/^(divr, F) is endowed with the norm 

l|p||-i/2,div,r •= IIPll-i/2,r + l|divrp||?-i/2,r- 

In H(curl,r2+), we can define the tangential boundary component that extends vr+u = 
1 / X u X 1 / and is bounded and surjective onto H~^/^(curlr, F), where curlp is the surface 
curl. Two traces, 7 “ and vr”, can also be defined from the interior domain. For functions 
u, V that are smooth enough on both sides of F with bounded support, we have the 
integration-by-parts formula 

(u, V X v)jr3\p - (V X u,v)]R 3 \r = ( 7 ^u, 7 r;:v) - ( 7 +u, 7 r+v), 

where (•, •)R 3 \r is the L^(M^ \ F) inner product and the angled bracket is the L^(F) in¬ 
ner product. This formula can be extended to u, v G H(curl,M^ \ F), with the angled 
bracket (•, •) now denoting the duality product of H“^/^(divr, F) and H“^/^(curlr, F). 
As is well understood since 0 El m, these two spaces are dual to each other with 
the duality product extending the inner product of square integrable tangential vec¬ 
tor fields. The rotation operator ^ ^ x acting on square integrable tangential 

vector helds can be extended to an isometric isomorphism between H“^/^(divr, F) and 
, F). Of capital importance will be the following jump operator, defined in 
H(curl2, M3 \ r) = {u G H(curl, M3 \ T) : V x V x u G L2(M3 \ T)}: 

[uJat := 7 “V XU- 7 +V X u. (2.7) 

Dependence with respect to time. Let us set up the kind of notation we will hence¬ 
forth use for functions of the space and time variables. We will assume that functions 
are defined as / : M —)■ X, where X is a Sobolev space on a domain or on F. The time 
derivative of / will be denoted /. We will also use the antidifferentiation symbol 

d-V(t):= f f{r)dr, 

Jo 

where we are using Bochner integration in the corresponding space X. A key space will 
be 


>V^(X) ;= {/ : M ^ X : / G C^-'(M; X), / ^ 0 in (- 00 , 0), G L\R- X)}, 

where the highest order differentiation is taken in the sense of X-valued distributions over 
M. In this space we can define the cummulative seminorms 

^ t 

Mf,t\X):=T [ \\f‘>(T)\\xdT, t>0. (2.8) 

£=0 

Note that / G W^(X) for k > 1, then is uniformly bounded. Another key 

space will be 

C^(X):={/gC^(M;X) : / ^ 0 in (-C50, 0)}. 
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The main theorems of this paper concern the dependence of the solution of fl2.6p 
and of = «Sc * with respect to (3 := —u x E‘“|r x i/, as well as the error of the 
approximations of to J and E^ to E. The theorems will only appeal to the fact that it 
is a closed subspace of H“^/^(divr, T) and will not break down even when Xh is inhnite 
dimensional. In particular the first theorem can also be read as a mapping estimate for 
the inversion of the continuous equation fl2.4p . which we rephrase as Corollary 12.31 The 
expression independent of h has to be understood as independent of the choice of Xh. 

Theorem 2.1. Let f3 G >V+(H“^/^(curlr, T)), let 3^ be the solution of 02.611 and E^ = 
Sc * J^. Then, e (L2(fi+)) n (H(curl, G (H-i/2(divr, T)), and there 

exists C > 0, independent of h and t, such that 

||J^(f)||_i/2,div,r + ||E^(^)||curi,o+ < C'max{c,c"^}iJ2(/3,t|H~^/^(curlr,r)) 

for all t > 0. 

Theorem 2.2. Let J and 3^ be the respective solutions of 02.4p and 02.611 and let E = 
Sc* 3 and E'^ = Sc* 3'^. Let Uf, : H-^/2(divr, T) ^ Xh be the orthogonal projection onto 
Xh- //J G W^(H“^/^(divr, T)), then there exists C > 0, independent of h and t, such 
that 

II J(t) - J"(t)||_i/2,div,r + ||E(t) - E'^(f)||,,ri,o+ 

< Cmax{c, c-2} H2{3 - Ilh3, t \ H~i/2(divr, T)) 


for allt>0. 

Corollary 2.3. Let /3 G W+(H~^/^(curlr, T)) and 3 be the unique solution ofVc*3 = 13. 
Then 3 G (H“^/^(divr, T)), «Sc * J G C^(L^(r2+)) fl (H(curl, 12+)), and there exists 
C > 0, independent oft, such that 

II J(^)l|-l/2,div,r + ||(5c * J)(t)||curi,o+ < C'max{c,c H2{f3,t\El ^/^(curlp, T)) 
for allt>0. 

Corollary 2.4. Let3 G >V2(H-i/2(divr,r)). ThenSc*3 G Ci(L2(f2+))nC° (H(curl, 12+)), 
and therefore Vc * J G C+(H~^/^(curlr, T))), and there exists C > 0, independent of t, 
such that 

ll(Ve*J)(t)|| -1/2, curl, r + II (‘5c * J)(^)||curi,o+ < CTnax{c,c ^}H2{3,t\H ^/2(divr,r)) 
for allt>0. 


3 Distributional potentials 

In this section we pave the way for the proofs of Theorems 12.11 and 12.21 by giving a very 
precise mathematical description of the problems: (a) given /3, compute E^ and J^, and 
(b) given J, compute E^ — E and — J. We will handle both problems simultaneously. 
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The theory is going to be developed with considerable latitude in the choice of the space 
Xfi. From now on, we will only assume that is a closed subspace of H“^/^(divr, F). 
The two limiting cases X^ = H“^/^(divr, F) and X^ = {0} will be discussed at the end 
of this section. We will use the polar set oi Xh 

XI := {r; G H-V2(curlr, F) : (^, ry) = 0 G X,}. 

Since Xh and are closed, we can dehne bounded operators : H~^/^(divr, F) 
H~^/^(divr, F) and Px° : H"^/^(curlr, P) —)■ H~^/^(curlr, P) such that 

Xh Px^^ = 0 and rj e XI = 0- (3.1) 

What these projections are is unimportant. The orthogonal projections on the orthogonal 
complement of the respective spaces can be taken for this role. 

The transient single layer potential for Maxwell’s equation can be rigorously dehned 
using a Laplace transform. The techniques are well-known and available in the literature. 
We will only introduce the definitions essential for a correct handling of the potential and 
its tangential component on the boundary. The following theorem lays the framework for 
causal vector-valued distributions, relating some distributions to their Laplace transforms. 

Theorem 3.1. [2ll Chapter 3] Let X be a Banach space and let f be an X-valued dis¬ 
tribution in M. The following statement on f 

there exists a continuous function g : M. ^ X such that g{t) = 0 for all t < 0 
and such that ||5'(t)|| < Cf^ for all t > 1 with m > 0, and there exists a 
non-negative integer k such that f = g^^'^ 

is equivalent to 

f admits a Laplace transform F = C{f} defined in C+ := {s G C : Res > 0} 
and satisfying ||F(s)|| < C'F(Res)|s|^ for all s G C+, where /i G M and C-p : 

(0, cxo) —)■ (0, cxo) is non-increasing and satisfies Cp{a) < Ca~^ for all a < 1 
for some C* > 0 and i >0. 

Following m, the set of all causal distributions characterized by Theorem 13 .1 1 will be 
denoted TD(X) (TD as in time-domain). Note that if / G TD(X) and A G B{X,Y), 
then Af G TD(r). Note also that for >V^(X) C TD(X) for all fc > 0. 

The Maxwell single layer potential. For brevity, we name the space 

M := H(curl, M^) n H(curl^ R^\T). 

We endow M with the norm 

I|u||m := I|u||k 3 + ||V X nllla + ||V x V x 

and note that ^ M —)■ H~^/^(curlr, F) and [-Jx : M —H~^/^(divr, F) are bounded. 
Given ^ G TD(H-P2(divr,F)), we can find a unique n G TD(M) such that 

c“^u -|-VxVxu = 0, [njx = —(3.2) 
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The first equation in fl3.2l) has to be understood in the sense of L^(M^ \ r)-valued dis¬ 
tributions, while the second one is an equation in the sense of H~^/^(divr, T)-valued 
distributions. Implicit to the fact that u is M-valued is the equality 7r“u = tt+u. The 
operator that given ^ outputs u is a convolution operator and we denote it as u = «Sc * 
We then denote Vc* ^ ■= tTt^c * which is a causal H“^/^(curlr, r)-valued distribu¬ 
tion. Note that the convolution product of causal distributions is always well defined, 
and therefore, the convolutions Sc * ^ and Vc * ^ can be extended to arbitrary causal 
H“i/^(divr, T)-valued distributions However, the relation of these operators to the 
transmission problem fl3.2p is unclear when ^ does not admit a Laplace transform. The 
existence and uniqueness of solution of fl3.2p is proved using the Laplace transform (see 
[1] for a recent and careful exposition). We now sketch the idea with wave speed c = 1. 
Given rj E H~^/^(divr, T) we look for U G M satisfying 

s^U + V X V X U = 0, [U]jv = -ST] (3.3) 

or, equivalently, we look for U G H(curl, M^) satisfying 

s^(U, v)]r3 (V X U, Vv)r3 = -s(? 7 , TT^v) Vv G H(curl, M^). (3.4) 

The solution operator for fl3.3p is a bounded linear operator S(s) : H~^/^(divr, T) ^ M. 
Proving a bound for ||S(s)|| in the style of the Laplace bounds of Theorem 13.11 we can 
show the existence of Sc G TD(H(H“^/^(divr, T), M)) such that S(-/c) = £{«Sc}. Then 
«Sc * ^ is characterized by its Laplace transform S(s/c)£{^}(s). 

A transmission problem. Let (3 G TD(H ^/^(curlr, T)) and ^ G TD(H ^/^(divr, T)). 
We look for u G TD(M) satisfying 

W^ii -|-VxVxu = 0, (3.5a) 

tTt-u — (3 E (3.5b) 

[uJat- f-G X/j. (3.5c) 

Let us first clarify where these equations take place. Equation fl3.5ap is an equality 
of L^(M^ \ r)-valued distributions. Equations fl3.5bp and fl3.5cp can be understood as 
equalities in the sense of distributions with values in H“^/^(curlr, T) and H“^/^(divr, T) 
respectively. In other words, we write the equivalent equations using the projectors fl3.ip : 

Px°(7rrU-/3) = 0, Px;,(HAf - c”^i) = 0. 

We remark that these are the equations for distributions of the time variable. An informal 
way of understanding it would be to assume that all quantities are functions of t and the 
equations are satisfied for all t. 

Proposition 3.2. For any f3 E TD(H"P2(curlr, P)) and ^ E TD(H"P2(divr, P)), prob¬ 
lem fl3.5p has a unique solution u G TD(M). 













Proof. The proof follows from a simple combination of arguments in [15] and [T]. We just 
sketch the main steps. We let B = C{j3} and S = For s G C+ we solve the coercive 

variational problem 

U(s) G H(curl, M.^), 7r,U(s) - B(s) G XI, (3.6a) 

(V X U(s), V X v)]R3 + (s/c)^(U(s),v)]R3 =-(s/c)(S(s),7r^v) Vv G Hq, (3.6b) 

where Hq := {v G H(curl, M^) : tt^-v G XI}. The solution of 03.61) . parametrized in the 
variable for the Laplace transform s, is the Laplace transform of an M-valued distribution 
u that solves 03.51) . Following techniques in [1] it is simple to prove that 

If |||U(s)||™,i,k 3 + ||V X V X U(5)||K3\r < C'(Re(f))|f|3(||S(s)||_i/2,div,r + ||B(s)||_i/ 2 ,curi,r), 

where C'(cr) = (^/(crminll, a^}) and C depends only on F. Theorem 13.11 can then be 
invoked to prove that U = /l{u} where u G TD(M) solves 03.5p . □ 

The solution of 03. 6 p can be written as U(s) = G/i(s)B(s) + E/i(s)S(s), using two 
bounded operators Gh{s) G ]y[^^ E/i(s) G ;B(H“^/^(divr, F), M). Tak¬ 

ing the inverse Laplace transform, this formula becomes the sum of two convolutions 

u = Gh* f3 + Sh* (3.7) 

where Gh G TD(B(H-F2(curlr, F), M)) and Sh e TD(;B(H“F2(divr, F), M)). By unique¬ 
ness of solution to 03.5p . if G TD(H“^/^(divr, F)) satisfies E Xh, then £h* = 0. 

Therefore 

^ * (^ — B/i^). (3.8) 

The next two results express the solution to the problems of Section [2] using the convolu¬ 
tion operators in 03.70 . In particular f3 ^ Gh * P will be the semidiscrete EFIE solution 
operator f Theorem 12. ip . while ^ Sh* ^ will be related to the associated error operator 
iTheorem 12.21) . 

Proposition 3.3. Let (3 G TD(H ^/^(curlp, F)) and define 

E^ = Gh*f3, j'' = -ca-'[E'^];v. 

Then G TD(H“^/^(divr, F)) and G TD(M) are characterized by 

J^EXh, Vc*J^-(3eXI E^ = S,*3} (3.9) 

Proof. This is a straightforward consequence of the definition of Gh and of the charac¬ 
terization of the Maxwell single layer potential by the transmission problem 03.2p . Note 
that the second equation in 03.91) can be equivalently written as a Galerkin semidiscrete 
equation: Vc * J^) = f3) for all G Xh- □ 

Proposition 3.4. Let J G TD(H define 

E^ = Sc*J-Sh*3, =-ca-^lE'^liv. 
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Then G TD(H ^/^(divp, F)) and G TD(M) are characterized by 


3’^ex^„ Vc*{J^-J)eXl E^ = S,*J\ (3.10) 

Therefore, ifE = Sc* 3, then 

E-E^ = Sf,*3 = SH*{3- n;,J), 3-3^ = -cd-^lE - (3.11) 

Proof. This result recasts Galerkin semidiscrete equations from the point of view of the 
exact solution: Vc * J^) = (/^^, Vc * J) for all G X^. Note that fl3.8p allows us to 

substract II/jJ in the argument of the operator of convolution by Sh- n 

Remark. The convolution Gh*(3+£h*^ is well defined for any causal pair of distributions 
(/3, ^) (by causal, we mean that their support is contained in [0, cxd)). This is due to the 
fact that the operator-valued distributions Qh and Eh are themselves causal. Moreover, 
Vi = Gh* (3 + Eh* ^ ^ solution of fl3.5p . because we can understand each of the equations 

in fl3.5p as the result of applying operators to the distributions Gh and Eh- What is not 
guaranteed is the uniqueness of solution of fl3.5p unless we restrict the space of possible 
solutions. We have opted for the set TD(M), which is somewhat restrictive but large 
enough for our purposes. Uniqueness can also be asserted in more general subspaces of 
the space of causal distributions, defined by the existence of a Laplace transform and 
some bounds on its behavior. 

Two particular cases. Assume that we take Xh = H“^/^(divr, F) so that X°^ = {0}. 
Then Eh = 0 and we are solving the problem 

Vc * J = /3, u = «Sc * J, 

i.e., we are dealing with the non-discretized inverse, also known as the continuous stability 
estimate. If we take Xh = {0} instead, then X)) = H~^/^(curlr, F). In this case Gh = 0, 
Eh = Sc, and we are dealing with 

u = 5c * I, 

i.e., we are handling the single layer potential. 


4 Proofs of the main results 

A first order system. Given (3 G TD(H-F2(curlr, F)) and | G TD(H-F2(divr, F)), 
we look for u G TD(M) and v G TD(H(curl; M^\F)) satisfying 

u — cV X V = 0, V + cV X u = 0, (4.1a) 

7r.u-/3GX°, [^^v]-|GX,. (4.1b) 

It is clear that if (u, v) solves fl4.ip . then u solves fl3.5p . 
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A strong version of the first order system. Most of this section will consist of the 
analysis of a problem related to (I4.ip . but written in M+ := [0,oo), using classical time 
derivatives and vanishing initial conditions. Our data are now smooth enough functions 
(3 : [0, cxo) —)■ H~^/^(curlr, F) and ^ ; [0, oo) H~^/^(divr, F), and look for 


u G C(M+;H(curl;M3)) nCi(R+;L2(M3\r)), 

V G C(R+;H(curl;R3\F)) nC^(R+;L2(R=^)), 

(4.2a) 

(4.2b) 

such that for all t > 0 


u(t) — cV X v(t) = 0, v(t) + cV X u(t) = 0, 

7r^u(t) - /3(t) G XI, [7rv(t)] - $,{t) G Xh- 

(4.2c) 

(4.2d) 

and 

u(0) = 0, v(0) = 0. 

(4.2e) 

An unbounded operator. Consider Pi ;= L^(R^\F) x L^(R^), equipped 
ural norm, and the spaces 

with its nat- 

U,, := {u G H(curl,R3) : vr^u G XI}, 

Xh ■■= {v G H(curl,R3\F) : [y.v] G Xh). 



In the domain D{A) := \Jh x ^h) we dehne the operator A ; D{A) Ti given by 
A(u, v) := (cV X V, —cV X u). 

Proposition 4.1. The operator A : D{A) -^TL is the infinitesimal generator of a unitary 
Co-group of operators in Ti. 

Proof. According to m Chapter 1, Theorem 4.3] or [TTl Theorem 4.5.4], we only need 
to prove that ±A are maximal dissipative. We first prove that 

(AW,1T)^ = 0 yWeD{A). 

For all W = (u, v) G U/i x V/j, we can write 

(AW, W)n = ((cV X V, -cV x n), (n, v))^ 

= (cV X v,n)R3\r - (cV x n,v)R3 = -c([ 7 ^v], vr^u) = 0, 

since [ 7 tv] G Xh and G 

Let X ; D{A) PL he the canonical inclusion of D{A) into Pi. We now prove that 
X — A ■. D{A) Pi is, surjective, i.e., given f G L^(R^\F),g G there exists 

(u, v) G D{A) satisfying 


u — cV X V = f, V + cV X u = g. (4.3) 

To prove this, we solve the following coercive variational problem 
u G Uft, (u, w)iR3 + c^(V X u, V X w)ffi3 = (f, w)ffi3+c(g, V X w)ffi3 Vw G (4.4) 
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and then define v = g — cV x u G We claim that (u, v) G D{A) and that fl4.3p 

is satisfied. To prove this, we first choose an arbitrary z G (T>(M^\r))^ C U/i as a test 
function in fl4.4p and show that 

(V X V, z)(x>'(R3\r))3x(n(R3\r))3 = (v, V x z)r3 

= (g - cV X u, V X z)r 3 = c"^(u - f, z)r3. 

This implies V x v = c“^(u — f) G L^(M^\r), and therefore v G H(curl, M^\r). For 
w G U/i, we have 

-([7rv],7r^w) = (V X v,w)R3\r - (v, V x w)r3 

= c“^(u - f, w)r 3 - (g - cV X u, V X w)r3 = 0. 

The observation nAJh = leads to 

(l7.vl,C) = 0 VC ex;. 

and therefore [ 7 tv] G Xh, which proves that v G V/j. The surjectivity of X — and 
dissipativity imply maximal dissipativity of A. 

To prove the X + is surjective, we solve a similar variational problem 

U G U/,, (u,w)r3 +C^(V X U, V X w)r 3 = (f, w)r 3 -c(g,V X w)r3 Vw G U/j 

and then define v = g + cV x u. The rest of the analysis is a minor variation of the 
previous case. □ 

Lifting operator. The next step is the construction of a lifting operator that will 
eliminate the non-homogeneous transmission conditions fl4.2dp . 

Proposition 4.2. Given /3 G H ^/^(curlp, T), ^ G H ^/^(divr,r), there exists a unique 
pair (u, v) G H(curl, R^) x H(curl,R^\r) satisfying 

u = cV XV, V = —cV X n, (4.5a) 

TTrU-lBeXl, Ijr^l-^eXh. (4.5b) 

Moreover, there exists Cy > 0, independent of the choice of X^, such that 

||u||curl,R3 + II v||curl,R3\r < Cp max{c, C ^}(||/3||-l/2,curl,r + 11^11-1/2,div,r) (4.6) 

Proof. Since Tir ■ H(curl, R^) —)■ H“^/^(curlr, T) is bounded and surjective, we can use a 
bounded right-inverse to build G H(curl, R^) such that 

TT^n^ = (3, ||u^|lcurl,R3 < Cl ||/3|| _i/2,curl,r- (4.7) 

Similarly, since [ 7 ,-] ; H(curl, R^\r) F) is bounded and onto, we can 

choose G H(curl, R^ \ F) satisfying 

[7rV^] = II V^||curl,R3\r < C 2 1|^ || _l/ 2 ,div,r- (4.8) 
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In the proof of Proposition 14.11 we have shown that X — ^ is surjective. Therefore, we 
can hnd (u°, v°) G D{A) satisfying 

u° — cV X v° = f := cV X — u^, (4.9a) 

v° + cV X u° = g ;= —cV X — v^. (4.9b) 

It is then obvious that (u, v) = (u^ + u°,v^ + v°) is a solution of fl4.5p . Moreover, 
proceeding as in the proof of Proposition 14.11 it is evident that u° is the unique solution 
of fl4.4p . Choosing w = u° as test function in fl4.4p and using fl4.7p - fl4.8p . we can bound 

min{l,c}||u°||,u,i,R3 < ^||u 0|123 + c2||V x 

< \/2max{l,c} maxICi, C2}(||/3|| -1/2, curl,r + 11^11-1/2, div,r)- 
Using now fl4.9p . we can bound 

||v°||curi,M3\r < ||v°||r3 + ||V X v°||R3\r < ||g - cV X u°|1r3 + ||c"^(u° - f)|lR3\r 

< \/2niax{c"\ l}^||u°|1^3 + c^HV x + c"^||f||R3\r + ||g||iR3 

< 2niax{c,c"^}niax{Ci,C2}(||/3||_i/2,curi,r+ ||^||-i/ 2 ,div,r) 

+ \/2niax{c, c"^}Ci||/3||_i/2,curi,r + \/ 2 C 2 ||||| — l/2,div,r* 

Hence, the estimate fl4.6p follows readily. Note that this estimate proves uniqueness of 
solution of fl4.5p . □ 

We are now ready to use the previous arguments to prove existence and uniqueness of 
solution to the non-homogeneous transmission problem 04.21) . We will need the spaces 

Cj(K+; X) := {/ e C‘(K+; X) : /«l(0) = 0, i<k-l}. 

and the following auxiliary result from the theory of evolutionary equations on Banach 
spaces (See [IHl Chapter 4, Corollary 2.5]). 

Theorem 4.3. Let A : D{A) G LL ^ LL be the infinitesimal generator of a contractive 
Co-semigroup of operators in LL and let F G Cq(M+;'H). The initial value problem 

U{t) = AU{t) + F{t) yt > 0, f/(0) = 0, 

has a unigue classical solution U G C^(M+;'H) nC(M+;Zl(M)) and we can hound 

\mt)\\n< f \\F{T)\\ndT, \\umn< Tll^lUdr, Vt > 0. 

Jo Jo 

Proposition 4.4. For all (3 G (H-i/2(curlr, P)) and $, G hV^ (H-V2(divr, P)), prob¬ 
lem 04.21) with data (/3|m+, £|]r+) has a unigue solution and for all t>0, 

||u(t)||curi,R3 < 4Crmax{c,c"^} {H2{/3,t\ H"^/2(curlr, P)) + H 2 {^,t\ H~^/^(divr, P))) , 

||v(f)||curi,R3\r < 4Crmax{c,c"^} (ip2(/3,t | H"^/^(curlr, P)) + H2{^A\ H~^'^^(divr, P))) , 

where Cr is the constant of Proposition \4F2 . 
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Proof. We first prove the result with slightly smoother data: f3 G Cq(M+; H ^/^(curlp, F)) 
and ^ e C2(M+; H"F2(divr, F)). Let L : H“F2(curlr, F) x H-F2(divr, F) ^ H(curl, M^) x 
H(curl, M^\F) be the lifting operator dehned by Proposition 021 We then dehne IFtc(^) = 
L{f3{t), ^{t)) for all t > 0, and consider the function F = Wtc — Wrc = L{/3 — /3, ^ 

The hypotheses on /3 and ^ imply that IFrc £ Co(®+J H(curl; M^) x H(curl; M^\F)) and 
therefore F G CQ(M+;?f). At this moment, we use Theorem 14.31 and Proposition 14.11 to 
dehne the unique solution Wq G fl C(R+; of 

Wo{t) = AWo{t) + F{t), t > 0, Wo(0) = 0. (4.10) 

It is then easy to show that (u, v) = IFrc + hFo is the unique solution of fl4.2p . (Uniqueness 
follows from uniqueness of solution to fl4.10p .) 

Using Theorem 14.31 we estimate 

\\Wo{t)\\n < f ||IUTc(r) - IUTc('r)|kdr 
Jo 

< Cr max{c, c~^} (iFi(/3, f; H“^/^(curlr, P)) + iFi(|, t; H'^/^(divr, P))) . 
Using the bounds for the lifting operator again, we next bound 
\\W{t)\\H<\\WTcmn + \\WomH 

< 2Urmax{c, c“^} (Tfi(/3, t; H“^/^(curlr, P)) + Tfi(^, t; H“^/^(divr, P))) . 
Analogously, we obtain 

\\wmn<\\WTcmH + \\Womn 

< 2Urmax{c,c-^} (i/i(/3,t; H-^/2(curlr, P)) + i/i(^, t; H-^/2(divr, P))) . 
The bounds in the statement of the Proposition follow then from the identities 

I|uWIILi,r 3 = ||u(t)||K3\r + I|c”^v(f)|l 23 , ||v(f)||2^^iR3\r = l|v(t)|lK3 + ||c"^fi(f)||2 3^r. 

To prove the result for the weaker hypotheses in the statement we just need to use a 
simple density argument. Alternatively, we can use a variant of Theorem 14.31 where 
F G C(M+;'H), F(0) = 0 and F is integrable. □ 


Extension operator. Consider a continuous function / : [0, cxo) —)■ X and 


{Efm := 



t > 0, 
t < 0, 


which is a causal A-valued distribution. If ||/(t)|| < C(1 +1™'), then Ef G TD(X). 

Proposition 4.5. Let (3 G >V2(H-F2(curIr,P)),^ G (H-F2(divr, P)) and let (u, v) 
be the unique solution of fl4.2p with data (/3|r+, ^|r+). Then (Uu, Ev) is the solution of 
04.1 p with (/3, as data. Therefore, Eu = Qh * + £h * 
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Proof. We first note that the hypotheses on f3 and ^ imply that 

H2{(3, 1 1 H-^/2^curlr, T)) + H2{^, t \ T)) < C(1 + t^) 

for some C > 0. It is now clear from the above and Proposition 14.41 that E\ G 
TD(H(curl; M^\r)). It follows from rom fl4.2ap that 

d-^u e C(M+;H(curP;M3\r)) nC^(M+;H(curl; M^)) n C2(M+; L2(M3\r)). 

The fact that there exists some other C > 0 such that 

||'9"^u(t)||M < [ ||u(r)||curi,R3 dr + |lc"V(t)||cu,i,R3 < C{1 + t^) 

Jo 

implies Ed~^n G TD(M) and therefore Eu G TD(M). Since u(0) = 0 and v(0) = 0, we 
have 


■^{Eu) = Ell = E{cV X v) = V X (cEv), 

LLL 

— {Ev) = Ev = —E{c'\/ X u) = —V X (cEu). 

(JjL 

The transmission boundary conditions can be interpreted in the distributional sense nat¬ 
urally. □ 

Proposition 4.6. Let f3 G >V+(H ^/^(curlr, T)), ^ G >V+(H ^/^(divr,r)) and u be the 
solution to fl3.5p .hc.. u = Qh*^+Lh*^. Thenu G (H(curl; M^))nC^(L^(M^\r)), (9“^Vx 
u G (H(curl; M^\r)), (9 “^[u]7v £ (H~^/^(divr, T)) and there exists a constant Cf > 0 

independent of the choice of Xh such that 

||u(t)||curl,K3 ||c9"V X u(t)||curl,K3\r + || C [u] at (t) || _i/2,div,r 

< Cf max{c, c-^} {H2{(3, t \ H"^/^(curlr, T)) + H 2 {^, t \ H-^/2(divr, T))) Vt > 0. 

Proof. Dehning v = —C(9“^V x u, (u, v) is the solution to fl4.ip . By Proposition 14.51 and 
uniqueness of solution of fl4.ip . (u|r^, v|r^) is the solution to fl4.2p with data (/3 |r+,^|r+) 
and we have estimates as in Proposition 14.41 The bound for x u follows from 

V = —cd~^V X u. The bound for 5“^[u]Ar results from the boundedness of the operator 
[7,] : H(curl, \ p) -> R-V^^curlr, P). □ 

Proofs of all the results of Section [21. By Proposition [T31 and in Theorem 12.II 
can be represented as EP = /3, = —c5“^[E^]. Their bounds and regularity are a 

consequence of Proposition 14.61 with ^ = 0. Corollary 12.31 is a special case of Theorem 12.II 
when = {0}. Similarly, the bounds for E — E^ = * J and J — = —c d~^ [E — E^] 

in Theorem 12.21 follows from Proposition 14.61 with (3 = 0 and Proposition 13.41 Corollary 
12.41 is a special case of Theorem 12.21 when Xh = {0}. 
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5 Comparison with existing results 

In this short section we compare our results with those of [1] (which develop the techniques 
in [25]) and [22] , 

The following two results are consequences of estimates in [T]. On the one hand, [U 
Theorem 4.4(a)] is a frequency domain result about V“^(also in [161 Lemma 2]). We 
convert it to a time domain result using m Theorem 7.1]. 

Theorem 5.1. If we solve Vc * J = /3 when f3 G W+(H~L2(^(.urlr, T)), then 
J G (H”^/^(divr, T)) and there is a constant C independent oft such that for arbitrary 
t > 0, 

||J(i)||-i/2,div,r < C' 7 ^max{l,t^} [ ||/3^‘^n^)ll-i/2,div,r dr. 
i- +t Jq 

The representation formula for the error of Galerkin semidiscretization [H Formula 
(43)], the discrete stability estimate [H Lemma 4.8] and the boundary integral operator 
bound [H Formula (32)] can be used to provide a time domain bound [H Formula (44)] 
without explicitly displaying the constant’s dependence on time. With the same esti¬ 
mates and [121 Theorem 7.1] we can prove the following theorem, which admits a better 
comparison with our results. 

Theorem 5.2. Assume hi- is a polyhedron. Choose Xh to be the lowest order Raviart- 
Thomas elements on a surface mesh over the polyhedron. Let J and be the solutions 
to 02.41) and 02.61) . If J G W® there is a constant C independent of 

t such that for all t > 0, 

||J(t) - J"(t)||_i/ 2 ,div,r < G-^max{l,t"} f ||(J - n;,J)('^)(r)||_i/ 2 ,div,r dr. 

t + r Jo 

Compared to Corollary 12.31 and Theorem 12.21 these estimates require higher temporal 
regularity of the input data. As mentioned in the introduction, there is a notable loss 
when going back and forth through the Laplace domain. Theorem 15.21 is written for 
a polyhedral domain, but the techniques in [1] can be easily extended to any Lipschitz 
domain and any discrete subspace dehned on it. The constant grows polynomially in time 
and henceforth the result does not rule out long time polynomial growth of the error even 
when the input data is compactly supported, which we show not to happen. 

The following result is a special case of m Theorem 3.1] when p = 1. It ruled out 
the long time growth given compactly supported input data. Nevertheless, understanding 
the underlying theory it relies on presents a daunting challenge for anyone who intends to 
extend the result to a more general setting. For example, the approach of gives few 
clues to a possible discrete error estimate. In comparison with Theorem 12.11 the estimate 
applies to scatterers only with smooth boundaries, while our results cover the scatterers 
with general Lipschitz boundary. Our result also tolerates input data with rougher spatial 
regularity, p) versus 

Theorem 5.3. Assume C- is a bounded domain with C°° boundary. Given 

(3 G n compactly supported in [0,T], let J be the solution to 
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Vc*3 = (3 and define E = S^*J, then E e C° (HH^^+)) n C|(L2(1]+)), J e C° (Hi/2(r)) 
and the following bound holds for all t >0, 

||E(t)||Hi(r 2 +) + IIJ(^)llHi/ 2 (r) < Ct sup f||/3(T)||Hi/2(r) + ll/3('r)||H3/2(r)) • 

rG[0,T] ^ '' 
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